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We present a measurement of the top-quark mass (m t ) in pp collisions at yfs = 1.96 TeV using 
ti events with two leptons (ee, e/i, or fifi) and accompanying jets in 4.3 fb _1 of data collected with 
the DO detector at the Fermilab Tevatron collider. We analyze the kinematically underconstrained 
dilepton events by integrating over their neutrino rapidity distributions. We reduce the dominant 
systematic uncertainties from the calibration of jet energy using a correction obtained from ti events 
with a final state of a single lepton plus jets. We also correct jets in simulated events to replicate the 
quark flavor dependence of the jet response in data. We measure mt = 173.7 ± 2.8 (stat) ± 1.5 (syst) 
GeV and combining with our analysis in 1 fb _1 of preceding data we measure rrij = 174.0±2.4 (stat)± 
1.4 (syst) GeV. Taking into account statistical and systematic correlations, a combination with the 
DO matrix element result from both data sets yields m t = 173.9 ± 1.9 (stat) ± 1.6 (syst) GeV. 

PACS numbers: 12.15.Ff, 14.65.Ha 



The masses of fundamental fermions in the standard 
model (SM) are generated through their interaction with 
a hypothesized scalar Higgs field with a strength given 
by a Yukawa coupling specific to each fermion species. 
The Yukawa coupling of the top quark corresponds to 
unity within uncertainties, and this value is constrained 
by a measurement of the top-quark mass (mt). In direct 
searches at the LHC for the standard model Higgs bo- 
son, both the CMS and ATLAS experiments observe local 
excesses above the background expectations for a Higgs 
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boson mass (m#) of approximately 125 GeV/c 2 [1, 2], 
decaying to diboson final states. Combined results in 
searches from the CDF and DO experiments at the Teva- 
tron show evidence for events above background expec- 
tation in bb final states [3]. It is therefore important 
to sharpen the measurement of m t , as its precise value, 
along with the mass of the W boson (mw), constrain the 
standard model prediction for ran through well defined 
radiative corrections. 

In pp collisions, top quarks (t) are primarily pro- 
duced in ti pairs, with each top quark decaying with 
BR(t -> Wb) ~ 100%. These events yield final states 
with either 0, 1, or 2 leptons from decays of the two W 
bosons. We consider the dilepton channels (2£) that con- 
tain either electrons or muons of large transverse momen- 
tum (px) and at least two jets. We analyzed such events 
previously [4, 5] using the neutrino- weighting (i/WT) ap- 
proach [6]. While the 21 channels have low background, 
the small decay branching ratio into leptons means that 



4 



m t measurements from these events remained statisti- 
cally limited unlike in channels with one lepton and four 
or more jets (£+jets). This situation has changed recently 
(e.g., Ref. [7]). Now, dominant systematic uncertainties 
from jet energy calibration, which have been larger [4] 
in the dilepton channel compared to ^+jets, are limiting 
precision of the m t measurement. In £+jets events, two 
quarks originate from W boson decay and yield a dijet 
mass signature that permits a precise calibration of jet 
energies for the measurement of m t in tt events [8] . While 
this calibration has greatly improved measurements in 
the £+jets channels, it has not been carried over to the 
calibration in other analyses. This is primarily due to 
differences in event topologies that can affect the details 
of the jet energy scale. 

We present a new measurement of m t using the DO de- 
tector with 4.3 fb _1 of pp collider data in the ee, e/i, and 
fifi final states. We improve the jet energy calibration for 
the accompanying jets using the energy scale from £+jets 
events [10]. Our approach differs from that of Ref. [11] 
in that we do not use the £+jets scale as a constraint in 
a combined fit of £+jets and dilepton events. Instead, we 
use this constraint as a calibration, and estimate the un- 
certainties of transferring that calibration to the dilepton 
event topology. This procedure demonstrates how the 
calibration obtained using the dijet constraint from mw 
can be applied to different final states, and has wide ap- 
plicability beyond the measurement of m t in 2£ events. 
We also employ flavor-dependent corrections to jet en- 
ergies for the first time in a dilepton analysis that sub- 
stantially reduce the uncertainties on jet energy resulting 
from jet flavor. The presented m t measurement is per- 
formed using the same data as Ref. [7], and is correlated 
with it as discussed below. 

The DO detector [12] is a multipurpose detector op- 
erated at the Fermilab Tevatron pp collider. The inner 
detector consists of coaxial cylinders and disks of silicon 
microstrips for track and vertex reconstruction. Eight 
layers of scintillating fibers arranged in doublets surround 
the silicon microstrip tracker and extend tracking mea- 
surements to forward pseudorapidities, r\ [13]. A 1.9 T 
solenoid produces a magnetic field for the tracking de- 
tectors. Uranium- liquid argon calorimeters surround the 
tracking volume and perform both electromagnetic and 
hadronic shower energy measurements. Thin scintilla- 
tion intercryostat detectors sample showers in the region 
between the central and end calorimeters. Three layers 
of proportional drift tubes and scintillation counters re- 
side outside the calorimetry, with 1.8 T toroids that pro- 
vide muon identification and independent measurement 
of muon momenta. 

We simulate tt events using Monte Carlo (MC) sam- 
ples for 140 GeV < m t < 200 GeV using the ALPGEN 
generator [14] and pythia [15] for parton fragmenta- 
tion. Backgrounds originate from Z/j* — > 2^+jets and 
WW/WZ/ZZ -> 2^+jcts production. For the former, we 
use ALPGEN combined with pythia, while diboson back- 
grounds are simulated entirely with pythia. We pass all 



MC events through a full detector simulation based on 
GEANT [16]. Backgrounds from instrumental effects that 
result in misidentified leptons are modeled using data. 

We use single and two-lepton triggers to select events 
for this analysis. Data and simulated events are recon- 
structed to provide the momenta of tracks, jets, and lep- 
ton candidates. Charged leptons are required to be iso- 
lated from other calorimeter energy deposits, and to have 
an associated track in the inner detector. Calorimeter 
and tracking information are combined to identify elec- 
trons. Track parameters in the muon and inner detector 
system are used to identify muons. We reconstruct jets 
with an iterative, midpoint cone algorithm with radius 
T^conc = 0.5 [17]. Jets are calibrated with the standard 
DO jet energy correction which is derived from data [18]. 
The method corrects the measured jet energy to the value 
obtained by applying the reconstruction cone algorithm 
to particles from jet fragmentation before they inter- 
act with the detector. We establish the efficacy of the 
method in the MC, where we compare the measured jet 
and the jet reconstructed from fragmentation particles. 
The jets in data and MC are calibrated independently so 
that their relative response is close to unity. This cor- 
rects for detector response, energy deposited outside of 
the jet cone, electronics noise, and pileup. The largest 
correction compensates for the detector response, and is 
extracted using 7+jet events in data and MC. We also 
correct jets for the px of any embedded muon and that of 
the associated neutrino. We initially apply this standard 
calibration [18] because it provides detailed px and rj de- 
pendent corrections. It also provides distinct corrections 
to jets and the imbalance in event transverse momentum 
($x) because several components (e.g., noise and out-of- 
cone effects) result from the jet reconstruction algorithm 
rather than any undetected energy. In the px range of 
jets found in tt events, the uncertainty of the standard 
DO jet energy calibration averages 2%, and is dominated 
by systematics. Because the flavor dependence of jet en- 
ergy calibration can yield one of the largest systematic 
uncertainties on our measurement [4], we have improved 
our analysis by accounting for this dependence. We use 
responses of single particles from data and MC to deter- 
mine the energy scale for different jet flavors. We correct 
MC jets by the ratio of data response to MC response 
according to their flavor to ensure that the MC reflects 
the flavor dependence in data, as in Ref. [10]. We calcu- 
late |?t as the negative of the vector sum of all transverse 
components of calorimeter cell energies and muon track 
momenta, corrected for the response to electrons and jets. 

Events are selected to have two leptons (ee, e^i, fj,fi) 
and two or more jets. The leptons must have px > 
15 GeV and the jets must have px > 20 GeV. Electrons 
and jets are required to satisfy \r/\ < 2.5, while muons 
must have < 2. We further require $x> 40 GeV in 
the channel. The efi events must satisfy Hx > 120 
GeV, where Hx is defined to be the sum of the pxs of 
jets and the leading lepton. In and ee events, we also 
require $x to be significantly larger than typical values 
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found in the distribution from Z boson events. These 
and all other selections are detailed in Ref. [19]. We ob- 
serve 50, 198, and 84 events with expected background 
yields of 10.4, 28.1, and 31.0 events in the ee, e/i, and 
channels, respectively. 

In £+jets events, one W boson decays to two quarks 
that fragment to jets. The invariant mass of this jet pair 
can be used to improve the calibration for all jets in these 
events. Complications arise because the four jets in the 
f+jets events can be incorrectly assigned to the initial 
four quarks. Energy from different partons is also mixed 
in the same jet due to a high jet multiplicity. Observed jet 
energies are also affected by color flow effects, which are 
different for the 6-quark jets and for jets from the decay 
of color singlet W bosons. These attributes are specific 
to a particular event topology such as £+jets. Neverthe- 
less, a scale factor based on the dijet invariant mass that 
is correlated with mw can be extracted. The most re- 
cent analysis of this kind by DO used 2.6 fb -1 of data and 
obtained a calibration factor of 1.013 ± 0.008 (stat) [10]. 
The uncertainty of 0.8% is smaller than that of the stan- 
dard jet energy correction and will decrease with addi- 
tional data. There are additional systematic effects on 
this energy scale that one must account for when ap- 
plying it to 6-quark jets in the £+jets analysis. These 
also affect our analysis, and we similarly evaluate the 
flavor dependence and residual energy scale systematic 
uncertainties directly on the measured m t to avoid dou- 
ble counting. These are quoted in Table II and discussed 
below. Beyond this, we have the possible difference be- 
tween 6-quark jets in dilepton events and &-quark jets in 
£+jets events and the effect of using a calibration based 
on a subset of the total data, each of which we discuss 
now in detail. 

The event topology is different in 21 and €+jets events. 
This has prevented significant progress in reducing the 
large standard jet energy scale uncertainties in dilepton 
analyses. To overcome this challenge and carry over the 
i?+jets calibration, we must account for the possibility 
that the energy scale of the 6-quark jets in the two chan- 
nels can differ. We calculate the energy scale, R 2e , for 
6-quark jets in the dilepton sample using responses for 
single particles that fall within the reconstructed jet cone. 
This is done by scaling single particle responses in MC 
to reproduce the energy response of jets in data [9] , giv- 
ing -R^ata' an d using particle responses from MC, giving 
i?jj C . We calculate the ratio of these two responses in the 
dilepton channel and the analogous ratio for 6-quark jets 
in the £+jets sample. The corresponding double ratio 

'Tjb ( b \ data 

2£(PT) "cr^)/<r s (^)' () 

varies between 1.001 and 1.003 depending on 6-quark jet 
Pt, Pt- The multiplicity of particles in 6-quark jets in 
£+jets events at the MC generator level is, after appli- 
cation of the offline jet algorithm, a few percent higher 
than in the dilepton sample, which is a sufficiently large 



difference to account for the observed value of IZ^. We 
therefore take 0.3%, the maximum excursion of TZ^ from 
unity, as a systematic uncertainty on carrying over the 
i?+jets scale to the jets in our dilepton sample. The £+jets 
scale is applied as a direct correction to the standard cal- 
ibration. 

The jet energy scale calibration obtained in Ref. [10] is 
based on a subset of the data, and we must therefore es- 
timate the effect of using the calibration on a larger data 
set. The instantaneous luminosity of the dilepton sam- 
ple is higher on average. We reweight the distribution of 
the number of primary vertices in the £+jets sample to 
match the distribution in the 4.3 fb _1 £+]ets data and 
recalculate the £+jets energy scale. This produces a neg- 
ligible effect. To account for a possible shift in the energy 
scale of the liquid argon calorimeter, we apply a correc- 
tion derived from 4.3 fb" 1 rather than 2.6 fb" 1 , and this 
yields a 0.7% shift in jet energy scale. From these stud- 
ies, we obtain a total uncertainty on the £+jets energy 
scale as applied to our analysis as the sum in quadrature 
of the statistical uncertainty (0.8%), 1Z\ t (0.3%), and the 
calorimeter calibration (0.7%). This yields a 1.1% uncer- 
tainty for applying the £+jets energy scale. 

The consequence of two neutrinos in dilepton events is 
an underconstrained kinematics. We employ the vWT 
technique to extract m t [6] due to its weak sensitivity to 
the modeling details of tt events. We integrate over the 
n distributions predicted for both neutrinos, solve the 
event kinematics, and calculate $t from the neutrino 
momentum solutions. The expected neutrino r\ distribu- 
tion in the dilepton channel is symmetric around 77=0 
and found to be well-described by a Gaussian distribu- 
tion. The width of the distribution decreases gradually 
with increasing m t (i.e., as the neutrinos become more 
central). Hence, we model the neutrino 77 distributions 
with a Gaussian probability distribution using a width 
parameterized as a linear function of m t . Several more 
sophisticated parametrizations were tested, but provided 
negligible improvement in expected precision in pseudo- 
experiments. By comparing the calculated $t to the 
measured $t for each event, we calculate a weight for a 
given choice of m t . For each neutrino rapidity sampling, 
we sum the weights calculated from all combinations of 
neutrino momentum solutions and jet assignments. We 
therefore arrive at a distribution of relative weight for a 
range of m t for each event. We found in Ref. [4] that 
most of the statistical sensitivity to mt is obtained from 
the first two moments of this weight distribution, the 
mean (fi w ) and RMS (<r w ). A coarse granularity of our 
sampling of the r\ distribution causes these moments to 
be unstable. To reduce this variation, we have increased 
the sampling for this integration by an order of magni- 
tude relative to our previous analysis [4]. This improves 
the expected statistical uncertainty on nit by 4%. Re- 
quiring the integral of this distribution to be nonzero 
excludes events with a measured $t that is incompatible 
with coming from neutrinos from tt decay. This intro- 
duces a small inefficiency for the tt signal and reduces 
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TABLE I. Parameters used to calibrate in the analysis 
of ee, e/t, and fifi channels and their combination. 



Channel 


Slope 


Offset [GeV] 


Pull width 


ee 


0.976 ±0.014 


0.03 ±0.16 


1.01 ±0.01 


e/i 


0.973 ±0.012 


0.43 ±0.14 


1.03 ±0.01 


nn 


1.038 ±0.022 


0.49 ±0.23 


1.06 ±0.03 



the background contamination in the final sample. Our 
final kinematically reconstructed data sample consists of 
49, 190, and 80 events in the ee, e/i, and /i/i channels, 
respectively. 

Probability distributions for \i w and a w are con- 
structed for background in each channel. Each back- 
ground component is normalized to its expected event 
yield. We generate distributions of ti signal probabil- 
ity as a function of ji w , a wi and m t . We use a binning 
that provides the minimum expected statistical uncer- 
tainty, as checked in pseudoexperiments. We perform 
a binned maximum likelihood fit to the probability dis- 
tributions, fixing the total signal and background yields 
expected in our data. The signal is normalized to the 
cross section calculated for ti production [20], evaluated 
at mt — 172.5 GeV. For all measurements, we obtain a 
likelihood (L) vs m t . We fit a parabola to the dependence 
of — In L vs. m t , and the fitted mass, mf 1 , is defined as 
the lowest point of the parabola. Point-to-point fluctua- 
tions mean that the initial placement of the window may 
result in an nonconvergent fit. We therefore iterate the 
fit around the current fit minimum. This results in a 
significant improvement in fitting efficiency, particularly 
in the dimuon channel. The final — In L vs. m t for data 
is shown in Fig. 1. The statistical uncertainty for each 
measurement is taken as the half-width of the parabola 
at 0.5 units in — InL above the minimum at m^'. 

The above procedure is followed for the extraction of 
m t from data and is used to calibrate the result as follows. 
We construct pseudoexperiments from signal and back- 
ground MC samples according to their expected yields 
and allow fluctuations in each such that the total equals 
the number of observed events. We perform 1000 pseudo- 
experiments for each channel, and measure in each. 
A linear fit of mf * vs the input nit provides a calibration 
for our method. We also calculate the pull width of the 
average estimated statistical uncertainty vs the rms of 
mf 1 values. The resulting slopes, offsets, and pull widths 
are given in Table I. The mf* and estimated statisti- 
cal uncertainty are corrected with these parameters. We 
obtain a calibrated mass measurement for the 4.3 fb -1 
sample in the ee, e/i, and /</i channels. 

The largest systematic uncertainties are associated 
with the jet calibration. We change the €+jets energy 
scale factor by ±1.1%, and perform our analysis to ob- 
tain a systematic uncertainty on m t of 0.9 GeV. The 
result of the £±jets analysis is a single scale factor aver- 
aged over all jet prs that are utilized in the dijet mass, 
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FIG. 1. -InL as a function of m t for the combined ee, e/i, 
and /i/t channels. A parabolic fit is shown near the minimum 
value in m t . 



i.e. dominated by light quark jets from W boson decay. 
As in Ref. [10], we estimate an uncertainty due to the dif- 
ference in pt distributions of 6-quark jets, in our case in 
dilepton events, vs the calibrating jets from the W — > jj 
sample. To estimate an uncertainty from this difference, 
we treat the px and n dependence of the uncertainty in 
the standard jet energy scale as a possible dependence 
of the residual energy scale following the calibration to 
£±jets. We calculate the average of the energy scale un- 
certainty for jets in the W — > jj sample. For each jet 
in the dilepton sample, we apply a shift corresponding 
to the difference between its uncertainty in energy scale 
and the W — > jj sample's average uncertainty in energy 
scale. Propagating this difference through the mass anal- 
ysis yields a 0.3 GeV uncertainty on m t . 

The flavor-dependent jet energy corrections described 
earlier provide MC-based mass templates that accurately 
reflect the data. As in [10], we propagate the uncertainty 
in these corrections and obtain a systematic uncertainty 
on m t of 0.5 GeV. The uncertainties due to flavor depen- 
dence and residual scale together with the uncertainty 
originating from the carry over of the jet energy scale 
from the £±jets sample account for the difference between 
6-quark jets in dilepton events and jets from W — > jj in 
£±jets events. 

We evaluate the effect of our uncertainty in model- 
ing initial state radiation (ISR) and final state radia- 
tion (FSR) by comparing two pythia samples having 
identical values of generated m t but different input pa- 
rameters taken from a CDF study [21] corresponding to 
an increased or decreased amount of ISR/FSR. Color 
reconnection uncertainties are estimated by comparing 
the analysis with PYTHIA Tune Apro and pythia Tune 
ACpro using [22]. Higher order QCD evolution is es- 
timated by comparing ALPGEN configured with pythia 
to mc@nlo with herwig [23] and this accounts for the 
uncertainty due to underlying event as well. To esti- 
mate sensitivity to uncertainties in the parton distribu- 
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TABLE II. Estimated systematic uncertainties on m t for the 
combined dilepton measurement in 4.3 fb _1 . 



Source 


Uncertainty (GcV) 


Jet energy calibration 




Overall scale 


0.9 


Flavor dependence 


0.5 


Residual scale 


0.3 


Signal modeling 




ISR/FSR 


0.4 


Color reconnection 


0.5 


Higher order effects 


0.6 


b quark fragmentation 


0.1 


PDF uncertainty 


0.5 


Object reconstruction 




Muon pt resolution 


0.2 


Electron energy scale 


0.2 


Muon p T scale 


0.2 


Jet resolution 


0.3 


Jet identification 


0.3 


Method 




Calibration 


0.1 


Template statistics 


0.5 


Signal fraction 


0.2 


Total systematic uncertainty 


1.5 



tion functions, we use CTEQ6M, and employ the method 
described in Rcf. [24]. 

We modify the jet energy resolution in MC events to 
reflect the resolution in data. We evaluate the effect of 
an uncertainty in this procedure on the extraction of nit 
by shifting the jet resolution by one standard deviation. 
We treat the electron and muon energy and momentum 
scales similarly and shift their calibrations within their 
uncertainties. 

Pseudoexperimcnts are used similarly to account for 
the uncertainty in the method that arises from the un- 
certainties on the offset and slope in the calibration of 
the fitted m t . We estimate the uncertainty due to the 
statistics employed in our templates of the ti probability 
distributions. We construct 1000 new templates, for both 
signal and background, and vary their bin contents within 



their Gaussian uncertainties. With these templates, we 
obtain 1000 new measurements from data and quote the 
rms of these values as a systematic uncertainty. We as- 
sign a systematic uncertainty on the signal fraction by 
shifting the background contributions in pseudoexperi- 
ments within their total uncertainty. 

We combine measurements in the three dilepton chan- 
nels using the method of "best linear unbiased estima- 
tor" [25]. We calculate each systematic uncertainty for 
the combined result, as given in Table II, according to its 
correlation among channels. The resulting measurement 
gives m t = 173.7 ± 2.8 (stat) ± 1.5 (syst) GeV. 

We combine this measurement with DO's measure- 
ment in the preceding 1 fb _1 of data using the i/WT 
and matrix weighting methods [4]. Some uncertainties 
evaluated in the 4.3 fb _1 sample are not available for 
the 1.0 fb _1 sample. We include the new uncertain- 
ties in the result from the previous analysis. We con- 
sider statistical uncertainties, as well as the following 
systematic uncertainties to be uncorrelated: calibration 
of method, template statistics, overall jet energy scale, 
and flavor dependence. We consider all other uncertain- 
ties to be fully correlated. The combined measurement 
yields m t = 174.0 ± 2.4 (stat) ± 1.4 (syst) GeV. This is 
consistent with measurements in other channels, and is 
the most precise single m t measurement in the dilepton 
channel to date. We have also improved the precision by 
combining the vWT results with the results of Ref. [7]. 
The statistical correlation of these two measurements is 
approximately 60%, calculated from pseudoexperimcnts. 
Accounting for this correlation, and correlations appro- 
priate to each source of systematic uncertainty, we obtain 
m t = 173.9 ± 1.9 (stat) ± 1.6 (syst) GeV. 
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